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We numerically determine the very rich phase diagram of mass-imbalanced binary mixtures of hardcore 
bosons (or equivalently - fermions, or hardcore-Bose/Fermi mixtures) loaded in one-dimensional optical lat- 
tices. Focusing on commensurate fillings away from half filling, we find a strong asymmetry between attractive 
and repulsive interactions. Attraction is found to always lead to pairing, associated with a spin gap, and to pair 
crystallization for very strong mass imbalance. In the repulsive case the two atomic components remain instead 
fully gapless over a large parameter range; only a very strong mass imbalance leads to the opening of a spin gap. 
The spin-gap phase is the precursor of a crystalline phase occurring for an even stronger mass imbalance. The 
fundamental asymmetry of the phase diagram is at odds with recent theoretical predictions, and can be tested 
directly via time-of-flight experiments on trapped cold atoms. 

PACS numbers: 37.1().Jk, 05.30.Jp, 71.10.Pm, ()3.75.Lm 



One-dimensional quantum liquids occupy a special place in 
the context of quantum many-body systems: indeed interac- 
tions of any strength lead to quantum fluctuations as strong 
as to discard Bose condensation for bosons and the Fermi liq- 
uid picture for fermions down to zero temperature. For suffi- 
ciently weak interactions a new unifying paradigm of the so- 
called Tomonaga-Luttinger liquids (TLL) emerges 1 1 1, char- 
acterized by the fact that all elementary excitations are gap- 
less, and both diagonal and off-diagonal correlations decay 
algebraically with the distance. Recent advances in the trap- 
ping of ultracold atoms in optical lattices allow to realize one- 
dimensional quantum liquids in a highly flexible way, with 
the possibility of fully controlling the statistics and the inter- 
action strength |2|. A series of recent experiments has demon- 
strated the physics of one-dimensional Bose gases with strong 
interactions up to the hardcore (or Tonks-Girardeau) limit |3|. 
A special role in the context of one-dimensional systems is 
played by binary mixtures, either bosonic, fermionic, or Bose- 
Fermi ones, for which TLL theory predicts the separation of 
spin and charge modes |[T1. In the case of particles with equal 
masses and repulsive short-range interactions, both charge and 
spin sectors can be gapless, and one recovers an effective pic- 
ture of two decoupled TLLs. Such a picture can be made un- 
stable via several mechanisms: via Mott localization in pres- 
ence of an underlying lattice and for integer total filling; via 
localization into a true (long-range-ordered) crystal (TC) in 
presence of a strong off-site repulsion; via phase separation; 
or via the formation of bound states {e.g. Cooper pairs for at- 
tractive interactions) leading to the appearance of a spin gap. 



Here we show that the TLL picture undergoes a complex 
series of instabilities in binary mixtures with mass imbalance 
between the two species. We focus here on quantum particles 
on a lattice with intraspecies hardcore repulsion and on-site 
interspecies interactions, describing at the same time spin- 1/2 
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FIG. 1 : Phase diagram of a mass-imbalanced atomic mixture with 
filling Ua = rib = 1/3. The phase boundaries are determined 
via QMC; the open diamonds indicate points in parameter space for 
which DMRG finds a fully gapped crystal phase TC2. The dashed 
line indicate the FK gap Ak/ Ja for kink-antikink pairs. The left 
panel presents a sketch of the phases: 2 TLL = 2 Tomonaga-Luttinger 
liquids; PS = phase separated; SG = spin gap; PSF = pair superfluid; 
CDW = algebraic charge density wave; TCI and TC2: true crystals. 

fermions, spin- 1/2 hardcore bosons, and mixtures of hardcore 
bosons and spinless fermions. The system Hamiltonian reads 



H 



E 



\^-{Jaa\ai+i + Jbb\bi+i +h.c.) + Uni^ant,b 



(1) 



in which a and h correspond to the two atomic species. Mass 
imbalance is controlled by the ratio j — J^/ J a and interac- 
tion by the ratio u — U j Ja- As far as the spectrum and the 
diagonal observables are concerned, we do not need to specify 
the statistics of the a and h operators; for what concerns off- 
diagonal observables, unless otherwise specified we will refer 
explicitly to hardcore bosons, satisfying bosonic commutation 
relations off-site and fermionic anticommutation relations on- 
site. The denomination of the many-body phases will also be 
mostly inspired by the case of hardcore -boson mixtures. 
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We consider both repulsive (u > 0) and attractive (u < 0) 
interactions, and we focus on the case of equal densities away 
from half filling, = rib = 1/p 7^ 1/2, p G N. The system 
with equal masses, j = 1, is integrable [41, and features spin- 
charge separation into two TLL in the repulsive case; in the 
attractive one, an a-b bound state appears, associated with the 
opening of a spin gap, and only the charge sector remains gap- 
less, giving rise to a paired superfluid (PSF) phase. The case 
of mass imbalance has been studied recently by bosonization 
||5}^71, and numerically via the density-matrix renormalization 
group (DMRG) and related approaches tSi- rlOJ . 

Making use of numerically exact methods, here we deter- 
mine comprehensively the rich phase diagram of the system 
with mass imbalance, j < 1, as shown in Fig.[T] The latter fig- 
ure refers to the case n — 1/3, but the qualitative features are 
generic for commensurate fillings. Our main findings are the 
following: 1) in the attractive case, the PSF phase is found to 
persist up to strong imbalance, at which the system becomes 
unstable to the formation of a TC of pairs. 2) in the repulsive 
case, on the other hand, the double TLL of the mass-balanced 
case survives up to a large imbalance, at which two different 
instabilities appear For sufficiently large repulsion, the two 
species phase-separate. For weaker repulsion, first a spin-gap 
opens, associated with fluctuating magnetic order captured by 
string correlations; this phase is the precursor of a TC phase 
with pinned b particles, and a particles localized in the inter- 
stitial regions. 

Before discussing the derivation of the phase diagram, we 
point out that the physics described in this paper is largely ac- 
cessible to current experimental setups on ultracold mixtures 
in one-dimensional optical lattices. The mass imbalance can 
be realized with heteronuclear mixtures [e.g. '"^K-^Li, for the 
fermionic case, ''^K-^^Rb for the bosonic case, and ^°K-*^Rb 
for the Bose-Fermi case 1 1 1 1), and with homonuclear mixtures 
in different hyperfine states, and it can be continuously tuned 
by using lasers with a wavelength close to the magic value or 
to an atomic resonance for one of the two species. The inter- 
species interaction can be tuned by Feshbach resonances, as 
widely demonstrated in the recent literature ifTTl . 

The phase diagram in Fig. [T] is the result of a joint nu- 
merical study based on quantum Monte Carlo (QMC) and 
DMRG|12|. Our QMC calculations are based on a canoni- 
cal formulation |T3l of the Stochastic Series Expansion ap- 
proach with directed loops [14\, applied to chains with up to 
L = 150 sites with periodic boundary conditions, and at tem- 
peratures /3Jb ^ L capturing the T = physics for both 
species. Our DMRG calculations apply to chains with up to 
i = 144 with open boundary conditions and retaining up to 
M = 1400 states. 

We start our discussion by the attractive case. For a 
very broad range of mass imbalance, the system displays a 
PSF phase, characterized by quasi-condensation of bound a- 
b pairs, giving rise to an algebraic decay of the pairing cor- 
relation function HI, Gab{r) = {a]b]bi^rai^r) ^ r^^/^f 



density-density correlations are also decaying algebraically as 

Kp I A cos(27rnr) . i 
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FIG. 2: Evidences for the absence of spin gap for moderate mass 
imbalance, (a) Green's function for both species (solid lines are fits to 
Aa(h-,d{r\L)-^/'^^'^'^^'')\ L = 96); (b) k„ coefficient (dashed lines 
are linear fits); all data are for u — 4; (c) Block entropy 5"; for chains 
with L — 120 sites (the curves for j — 0.5 and j = 0.7 are shifted 
to improve readability); thick lines are a reference for c = 2 (black) 
and c = 1 conformal field theories (red); (d) Fermionic spin-spin 
correlation functions, calculated via DMRG between the sites L/3 
and 2L/3 of a chain with L — 96. The data for j — 0.4 and j = 
0.2 have been reseated to improve readability (by a factor 10^^ and 
10^*, respectively). 

n-i = rii^a + ni,b, and Kp is the charge TLL parameter For 
equal masses, Kp > 1 1 1 1 for all u < 0, so that the dominant 
correlations are the pairing ones. Moreover A'^ — due to 
the presence of the spin gap (see below). We extract the Lut- 
tinger exponent from the slope of the density structure factor 



at 9 0, Sp{q) = J2r^^P('^1'^)'^P 



Kp qj-K, and we 



Cp{r) {riini+r) ~ {ni) {ui+r) ~ 



i here. 



find that mass imbalance leads to a reduction of Kp, consis- 
tently with what observed for other fillings in Ref. ifTOl . For 
large mass imbalance, Kp becomes smaller than one: this cor- 
responds to the loss of quasi-condensation, in favor of a quasi- 
solid phase (or charge density wave, CDW), with dominant 
density correlations. This phase is the precursor of a quantum 
phase transition to a TC of pairs - phase TC2 of Fig.[T]- with 
the onset of long-range density order at wavevector Q = 2?™. 
We determine the extent of the TC phase via QMC by deter- 
mining the mass imbalance at which Sp{Q) starts diverging 
linearly with system size; and by DMRG detecting the onset 
of the exponential decrease of Gab and Cp, marking the open- 
ing of a charge gap ifTSl . The TC instability is well under- 
stood coming from the Falicov-Kimball (FK) limit, Jf, = 0. 
In this limit, which reduces to a ID lattice gas in a static po- 
tential, we find that the ground state corresponds to the TC 
of pairs for all values m < 0. The gap Ak to the formation 
of kink-antikink pairs in the TC is found to be a non-monotic 
function of displaying an intermediate maximum. We ob- 
serve that, for small the boundary of the crystalline region 
follows closely the locus at which the gap equals Jb, and 
it has a re-entrant shape mimicking the non-monotonic be- 
havior of the gap as a function of u. This suggests that the 
quantum melting transition corresponds to a condensation of 
kink-antikink pairs in the ground state. 

The repulsive side of the phase diagram is more com- 
plex. Combining bosonization with renormalization group 
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calculations up to two loops, Refs. |5] |6] conclude that a spin 
gap should open for any infinitesimal mass imbalance as 
As w Aexp(-yl'/| Ja - M) for Jb < Ja, where A 
Ja, Jb- On the contrary, Ref. 7 also based upon bosoniza- 
tion, concludes that the spin gap is absent in the repulsive 
case. All our numerical findings point toward the persistence 
of a fully gapless TLL behavior for both the charge and spin 
sector over a dominant portion of the repulsive phase dia- 
gram. Our conclusion is based on a number of crossed evi- 
dences. First of all, we observe that the one body correlation 
functions Ga{r) = (a-ai+r) and Gb{r) = {b\hi+r) can be 
very well fitted with the simple power-law form Ga{b) [f) — 
A„(b)(i(r|i)-i/(2K„(„)^ where d{r\L) = L| sin(7rr/L)|/7r 
is the conformal distance (see Fig. |2ja)). In particular we 
find that, for weak and moderate repulsions, Ka_(b) > 0.5, 
which implies that the momentum distribution ng^(b){q) — 
exp(i(3'r)G'Q(f,)(7') displays a quasi-condensation diver- 
gent peak at g = 0, to be detected in time-of-flight experi- 
ments (see below). Moreover we exploit the fact that an ex- 
plicit counting of the number of gapless degrees of freedom 
in the system comes from the central charge c of the confor- 
mal field theory corresponding to our model of interest. This 
quantity can be directly extracted via DMRG, using the funda- 
mental result that the entanglement entropy (EE) of a bound- 
ary block of the system grows with the size I of the block llT6l 
as Si = -Tr(pi log2 pi) « f log^iAL)] + + const, 

where pi is the reduced density matrix of the boundary block, 
and open boundary conditions are employed. Fig. |2|b) (see 
also Ref. 1, 15i ) shows that the scaling of the EE is fully con- 
sistent with c = 2 in the repulsive regime, providing further 
evidence for the fact that the TLL has two gapless components 
even for a significant mass imbalance. 

Finally, using QMC we gain further insight into the gap- 
less phase by investigating the spin-spin correlation func- 



B 



cos(27rnr) , 



where 



tionCr^(r) = (5f5f+,) . 
5,f = (rij a — ni^b), and is the spin TLL parameter; in 
absence of a spin gap, > 1. We extract the Luttinger ex- 
ponent from the low-q behavior of the spin structure factor, 
So-{q) ~ q/n, giving the finite-size estimate Ka-{L) ~ 
[{L/2)Sa[2TT/L-L) + (i/4)5^(47r/L;L)]/2. At the SU(2) 
invariant point j — 1, Ka{L) is known to obey the Kosterlitz- 
Thouless (KT) critical scaling Ka{L) 1 + [a log(L/Lo)]~^ 
(Qj, which implies that the quantity K,a{L) ~ [K„{L) — 
scales linearly with log(L). Quite remarkably, this scaling 
law remains valid over a broad range of j < 1 values, and 
even for strong inter-species coupling, as shown in Fig. |2|c). 
This means that, even if the SU(2) symmetry is broken in 
the Hamiltonian, it appears to be restored in the ground-state 
of the system, namely that a moderate mass imbalance is 
an irrelevant perturbation and that, under a renormalization 
group treatment, the system flows to the SU(2)-symmetric 
fixed point Ka{L — > oo) 1 over a broad range of j val- 
ues. This picture is further corroborated by comparing the Ion- 
gitudinal spin-spin correlations Ca (r) and the transverse 
ones in the fermionic case, ci^^\r) — {SfSf_f^j.), where 
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FIG. 3: (a) Breakdown of SU(2)-invariant scaling of k^, for u = 2 
and decreasing j; dashed lines are fits to ai/ log(L) + log{L)'^^; 

(b) 6-hole density correlator for nearest neighbors, L = 30; the 
shaded areas mark the transition regions among the various phases; 

(c) String structure factor as a function of system size (same param- 
eters as in (a)); (d) Decay exponent of the string correlation function 
as function of j. 

Sf = albi+blaf, the two correlation functions show the same 
algebraic decay over a broad range of j values , and they even 
appear to coincide at short range, as shown in Fig. |2|d). 

On the other hand, a careful study of the behavior of the 
system shows that the two-component TLL becomes indeed 
unstable to the formation of a spin gap (SG) for strong mass 
imbalance (j < 1/3 for n — 1/3). The opening of a spin gap 
leaves a clear signature in the Luttinger exponent K^, which 
stops flowing to the SU(2) symmetric value: this is seen in 
the fact that k„ begins to scale /oifer than linearly in log(L), 
which is a manifestation of the fact that Ka-{L) flows to val- 
ues smaller than 1, as shown in Fig.|3ja). We use this criterion 
to locate the opening of the spin gap in Fig.[T] Moreover the 
breaking of SU(2) symmetric scaling is also seen in the spin- 
spin coiTelation functions, which acquire a different behavior 
at short and long distances (see Fig. |2]^d)). The opening of a 
spin gap signals the formation of a bound state of some sort: 
but what could this bound state be in a purely repulsive sys- 
tem? We argue that a bound state can appear between the a 
particles and the holes of the b species, which are in a com- 
mensurate proportion m = (1 — n)/n to the a particles. As 
shown in Ref. [17] a sufficiently strong mass imbalance can 
indeed bind together particle composites in ID quantum flu- 
ids - in the specific case of n = 1 /3 we have tiimers formed 
by one a particle and two b holes (6). For strong enough re- 
pulsion (u > 1) the apparition of a-2b composites is seen 
in the nearest-neighbor density correlations of the b holes, 
C„„j, — ((1 — ni,b){l ~ Ui+i^b)), which is enhanced for in- 
creasing mass-imbalance, as seen in Fig. |3jb). 

The a-2b composites have mutual repulsive interactions, as 
observed in Ref 18 so that their binding favors configura- 
tions of the type ..aObaObOaObOab... in which the a and b par- 
ticles alternate spatially. Such configurations, albeit lacking 
long-range density and magnetic order, contain strong non- 
local coiTelations that can be captured using a string correla- 
tion function^ 0,{r) = -(si e^^^t=^^l^l Sf+^) . As 
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FIG. 4: Evolution of the momentum distribution of both species 
on a L = 60 chain, for interactions spanning the repulsive and 
the attractive side at fixed mass imbalance. From top to bottom: 
j = 2/3, 1/10, 1/20. 

shown in Fig. [3];c-d), string correlations are significantly en- 
hanced when entering the SG phase: indeed the string struc- 
ture factor, Oz,int — X]r'^z('")' shows a more pronounced 
divergence with system size, which signals the slower decay 
of Oz{r). As shown in Ref. T?" in presence of spin-charge 
separation the dominant decay of the string correlations is 
governed solely by the Kp exponent, Oz{r) ^ r^^p; there- 
fore Kp controls the divergence of Oz.int ^ L^^^p, and in 
Fig. [3jd) we observe that the 1 — Kp exponent (extracted 
from the low-g behavior of Sp{q)) is significantly enhanced 
in the SG phase. Given that both the spin and charge struc- 
ture factor exhibit a peak whose divergence with system size 
is hindered by the slow size-scaling of the divergence expo- 
nent 1 — Kg ~ Kp (15], we find that the string correlations 
best characterize the "fluid" magnetic order of the SG phase. 

The bound trimers appearing in the SG phase undergo crys- 
tallization - phase TCI in Fig.[T|- for moderate repulsion and 
extreme mass imbalance. The TCI phase is analogous to the 
one discussed in Ref. [TS] (see also Ref. ITST ). This phase is 
marked by a suppression of C„„ which attains its absolute 
minimum C„„ ^ = 1 — 2n for j ^ 0, due to the suppression 
of the weight of configurations with contiguous trimers. The 
resulting non-monotonic behavior of C„„ as a function of j, 
shown in Fig. |3jb) (see also Ref. I.15J ), provides further ev- 
idence that the SG phase is a liquid of pre-formed trimers - 
the enhancement of C„„ f, in that phase is due to trimer bind- 
ing, while the suppression in the TCI phase is due to crystal 
ordering of the trimers. 

Our theoretical results have immediate consequences for 
current experiments on one-dimensional mixtures of mass- 



imbalanced cold atoms. Such experiments can probe both 
attractive and repulsive interactions, within the same exper- 
imental conditions, via the use of Feshbach resonances fTl]. 
Our phase diagram reveals a fundamental asymmetry between 
the attractive and the repulsive case for weak and moderate 
mass imbalance, with the formation of bound pairs on the 
attractive side and the absence of spin gap on the repulsive 
one. This asymmetry is very well seen in time of flight ex- 
periments probing the momentum distributions na{k), iii,{k), 
which are very broad on the attractive side, while they ex- 
hibit sharp quasi-condensation peaks on the repulsive side - 
as shown in Fig. |4] This asymmetry can be used as strong ev- 
idence of pairing on the attractive side. A partial symmetry is 
recovered only for strong mass imbalance, with the opening 
of a spin gap in the repulsive case, and the occurrence of a 
crystalline phase for both signs of the interaction. This is also 
well captured by the momentum distributions, showing this 
time a suppression of the quasi-condensation peaks on the re- 
pulsive side due to the appearance of the SG phase and of the 
crystalline phase. Finally the availability of high-resolution 
in-situ imaging ll20l sensitive to the spin 11211 allows to detect 
the formation of bound trimers and the enhancement of string 
correlations characterizing the SG phase. 
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SUPPLEMENTARY MATERIAL 

Attractive case: PSF-CDW and CDW-crystal transition 

In this section we focus on the attractive case u < 0. 
The transition from the PSF phase to the CDW phase, ob- 
tained for decreasing hopping ratio j, is characterized by 
an inversion of the hierarchy between the pairing correla- 
tions, Gabir) ^ r^^/^p^ and the density-density correlations 
Cp{r) ^ cos(27rnr) r^^f . In the PSF phase Kp > 1, while 
in the CDW phase Kp < 1. The Kp exponent, extracted from 
the low-q behavior of the density structure factor as discussed 
in the text, is shown in Fig. [5] as a function of j for various 
system sizes. A finite-size extrapolation allows to identify 
the critical j value at which Kp traverses the critical value 
Kp = 1. The composite a-b pairs form therefore a TLL 
with effective exponent KcS = Kp/2. In the CDW phase, 
we obtain a TLL with Kcs < 1/2, which is quite remark- 
able in a system with on-site interactions only. Indeed in one- 
component TLLs such a situation is realized only in presence 
of non-local interactions, such as nearest-neighbor or dipo- 
lar ones (R. Citro et ai. New J. Phys. 10 , 045011 (2008); 
T. Roscilde and M. Boninsegni, New J. Phys. 12, 033032 
(2010); M. Dalmonte et ai, Phys Rev Lett. 105, 140401 
(2010)) or when considering highly-excited states such as a 
super-Tonks-Girardeau gas (G. E. Astrakharchik et ai, Phys. 
Rev. Lett. 95, 190407 (2005)). In our case, although the in- 
teractions between the a and b particles are on-site only, the 
bound a-b pairs have effective longer-range interactions, de- 
caying exponentially with the distance. 

These effective interactions are responsible for the transi- 
tion to crystalline order (phase TC2) at extreme mass imbal- 
ance. This transition occurs at the critical value Kj^g — 
where p ^ 1/n G NUJ. Since in our case p — 3, one obtains 
Kp"^^ — 4/9 at the transition line, and Kp = in the crys- 
talline phase, marking the onset of long-range density-density 
correlations. This implies the appearance of a peak diverging 
linearly in L at wavevector Q — 27rn in the structure factor for 
the total density Sp{q), or alternatively in the density structure 
factor of both species 




0.5 -^r CDW 



TC2 



• 


L=48 


• 


L=60 


• 


L^72 


o 


L=96 


• 


L=150 



0.2 



0.4 



0.6 



0.8 



FIG. 5: Luttinger exponent Kp as a function of mass imbalance for 
the attractive case u = ~1 and n = 1/3. The shaded areas mark the 
estimated transition regions between the various phases. 
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FIG. 6: Density structure factor for a and 6 particles at = — 1 and 
various j values. Solid lines are fits to the form ai + 02/!/"^. 
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(2) 

The onset of a linear-in-L divergence of the structure factor 
peak is shown in Fig.|6] where the peak value divided by L is 
seen to converge to a finite value in the thermodynamic limit. 
This criterion has been used to determine the phase boundary 
between TC2 and CDW in Fig.[T] 

A further clear signature of the TC2 phase is the exponen- 
tial suppression of the pairing correlator Gab, in sharp contrast 
with the algebraic decay of both CDW and PSF phases. We 
have used DMRG calculations with up to L = 144 sites for 
fermionic a and h species in order to determine the extent of 
the TC2 phase (diamonds in Fig.[T]). A typical set of data for 



j = 0.03 is presented in Fig. |7] there it is seen that for it = — 1 
and —2 the correlator Gab decays exponentially; interestingly, 
for stronger attraction, i.e. for more tightly bound pairs, the 
correlation function acquires again an algebraic decay, show- 
ing the re-entrant nature of the TC2 phase. We observe that, 
close to the onset of the TC2, the decay of the Gab correlator 



is indeed consistent with an exponent K 



4/9. 
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FIG. 7: Pairing correlation function on a L = 144 system starting 
from I//6 for j = 0.03 and different interaction strength. 

Repulsive case: entanglement entropy analysis 

As mentioned in the main text, the entanglement entropy 
(EE) allows to extract the central charge of the system| 16 1. In 
Fig. |8] we consider two different approaches to extract c. In 
the left panel, we plot S{1) as a function of the block length 
I at a fixed hopping asymmetry j = 0.5 and fixed length L, 
for different values of u; in the repulsive case, the numeri- 
cal results are compatible with a completely gapless system, 
whereas in the attractive case the emergence of the spin gap 
reduces the central charge to 1 . In the right panel, we plot the 
EE of the half chain 

S{L/2) = ^ log2(L/7r) + 0{l/L) + const. (3) 

which, for L — 6m, m e N, does not exhibit any oscillatory 
correction. For u > 0, all datas fit very well with the CFT 
prediction with c ~ 2, whereas a clear bending at large system 
sizes indicates that a gap opens in the attractive case. 




FIG. 8: Entanglement entropy of a subblock of length I. Left panel: 
EE scaling for L = 120, j — 0.5 and different values of u; for 
u > and u < 0, results are compatible with c — 2 (black thick 
line) and c = 1 (dashed line) respectively. Right panel: scaling of 
EE for I = L/2 at different system sizes. In the repulsive case, datas 
are compatible with c — 2, whereas in the attractive one datas for 
larger system deviates toward c = 1 at large L; dashed lines are 
linear fits with c — 2.01, 2.07, 2.05 from top to bottom, whereas the 
black, thick line represents c — 2. Data points for (j — 0.5, u = 
1) and {j = 0.9, it = 2) have been shifted by 0.2, 0.5 to improve 
readability. 



Repulsive case: nearest-neighbor density correlations 



Fig.|9]shows the nearest-neighbor correlation function for h 
holes, C„„ ^ = ((1 — ni i){\ — nj+i {,)) in the repulsive case. 
The calculation is performed via exact diagonalization on a 
chain with i = 12 sites, which allows to access properly the 
true ground state in the limit j — >■ 0, both in the case of crys- 
talline order and of phase separation. The latter limit is hard 
to reproduce appropriately with QMC or DMRG simulations, 
due to the abundance of metastable states close in energy to 
the ground state. We see that C„„ 5 is very sensitive to the 
succession of phases traversed by the system. As explained 
in the text, for small u values C„„ is suppressed as j — > 
due to the formation of the crystal of a — 26 trimers, imply- 
ing that C„„ I has to attain its minimum value 1 — 2ri = 1/3. 
As j is reduced, two competing effects result in a non-trivial 
evolution of C„„ ^. The formation of a — 26 trimers leads to 
an enhancement of the C„„ 5 correlator: the kinetic repulsion 
between the h holes, which would be normally present in a 
ID single-component gas with hardcore repulsion, is strongly 
suppressed due to the decrease of the hopping ratio j', and 
therefore it is overcome by the binding effect of the light 
a particles, which gain energy by delocalizing over two (or 
more) 6 holes. Yet the trimers that form via this mechanism 
experience an effective mutual repulsion at sufficiently small 
u values, due to the kinetic repulsion between light a particles. 
This leads to the suppression of configurations with adjacent 
trimers; a complete suppression of such configurations would 
lead the system to the trimer crystal state, which, as mentioned 
before and in the text, coiTesponds to an absolute minimum of 
C„„ i- Given that trimers are only loosely bound for small u, 
the effect of the repulsion between trimers dominates over the 
effect of trimer formation, so that C„„ ^ is seen to decrease 
monotonically for decreasing j {u — \ data in Fig. [9]l. On 
the other hand, for larger w, the enhancement of C„„ 5 due to 
trimer formation dominates over the suppression due to trimer 
repulsion, so that C„„ is seen to increase for decreasing j 
down to the transition point to the crystalline phase, at which 
it drops down to the j — ^ limit (u — 2 data in Fig.|9]l. 

On the other hand, for even larger u the effective attrac- 
tion between h holes mediated by the a particles begins to 
involve more than two holes at once, as the a-h interaction 
overcomes progressively the kinetic repulsion between a par- 
ticles and favors their aggregation over extended clusters of 
holes. This phenomenon culminates into the phase separation 
between hole-rich and hole-free regions 1 10 18| upon decreas- 
ing j , namely upon suppression of the kinetic repulsion be- 
tween holes. This succession of phenomena leads to a much 
stronger enhancement of C„„ 5 as j decreases, with a steep 
increase when phase separation is achieved. When full phase 
separation is realized in the limit j 0, C„„ 5 takes the value 
[(1 - n)L - 1]/L (namely 7/12 = 0.583333 for L = 12, as 
seen in Fig.|9]l. 
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FIG. 9: Nearest-neighbor correlation function for b holes (b) for a 
chain or L = 12 sites. 



FIG. 10: Scaling of the peak in the spin structure factor for the repul- 
sive case, u = 2. 



Decay of string correlations from bosonization 

In the following, we will study the dependence on the TLL 
parameters of the generaUzed string correlation function: 



Ol{r) = St exp 



where we have defined Sf 



E 



Su 



S: 



i-\-r 



(4) 



2 . After applying stan- 
dard bosonization identities to a, b operators by introducing 
the bosonic 0^, (/if, fields| 1 1, one can define effective spin and 
charge fields as (j)cr,p — {(pa T 4'b) I such that: 



dr(j)a{r) 



-lV^^{r)V^^{T) 



(5) 



where V^p{r) — e^''-°"^'''p'^'^^ are vertex operators related to 
the charge and spin fields and 7 is a constant; additional con- 
tributions with higher scaling dimension play no significant 
role in the remaining, and can thus be neglected. The non- 
local contribution in Eq. |4]can be reduced to: 



exp[*7r/3 ^fc'] 

2<fc<i+r 



exp[z7r/: 



e■x.Y>[i^/^T(i[(j)a{^ + r) - cj)a{i)]]- 



In addition to a pure spin contribution which stems from the 
first term in Eq. (|5]l, and which decays faster than r^^, the 
string correlator contains a more complicated, second contri- 
bution, stemming from vertex operators, and which we will 
indicate as O^'^. Assuming spin-charge separation, O^'^ can 
be written as: 

Of{r) cx {V;^{l)V-^{^ + r))-K 

X {V^^{^)V-"^{i)V,P^{^ + r)V^^ {i + r) 
^ Mi^\r) * (l/r^")- 

whose asymptotic decay is affected by both charge and spin 
contributions: the charge part has a fixed exponent, whereas 
the spin one M''a\r) has a /3-dependent form; in particular. 



for P — 2 coiTelation considered in this work, the slowest 
decaying part of M''a\r) is a constant, so that the dominant 
term in the string correlation functions decays as: 



Ol{r) cx l/r 



(6) 



It is worth noticing that, from Eq. (|6]), one can easily recover 
standard spin correlations for = decaying as i/r^'^p+i^'' , 
The fact that the decay of string correlations is uniquely 
governed by the Kp exponent is indeed remarkable; as dis- 
cussed in the main text, this leads to a stiing structure factor 
diverging with system size as Oz,int ^ L^~^>' (and Kp < 1 
for all j values). This result has to be compared with the scal- 
ing of the peaks of the charge and spin structure factor, namely 
for Sp{Q) and Sa{Q) with Q = 27rn. These peaks come 
from the dominant oscillating term (going like cos{Qr)) of 
the Cp{r) and Ca{r) correlation functions, whose amplitude 
decays like r^Kp-K^ ^ TYiis means that in general the Q-peak 
height scales with system size as L^^^p^^" . In the SG phase 
the Ka exponent vanishes for L — > 00, so that, strictly speak- 
ing, the true asymptotic scaling of Sp{Q) and Sc{Q) is in- 
deed the same as that of Oz,int- But in practice the scaling of 
K„ (L) to zero is extremely slow - typically logarithmic, a fact 
which generates a huge finite-size correction to the effective 
scaling of Sp{Q) and Sc{Q). In our simulations we observe 
that these quantities grow more slowly than log(-Z^) over the 
entire SG phase for all the system sizes we considered (up to 
L = 96 - see Fig. 10 1. To see them diverge at all with sys- 



tem size one might need exceedingly big system sizes, both 
for simulations and for experiments. Therefore stiing coiTela- 
tions stand as a unique tool to characterize the correlations in 
the SG phase. 



Repulsive case: onset of the crystalline phase 

In this section we provide further details of the estimate of 
the transition line between the SG and the TCI phase in the 
repulsive case. The transition is obtained by a joint finite-size 
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FIG. 11: Upper panel; scaling of the superfluid density of the a par- 
ticles for u = 1.5. Lower panel: scaling of the peak in the density 
structure factor for the a particles. All solid lines are fits to the form 



scaling analysis. We focus on the scaling of the peak in the 
structure factors Sa{b) io) appearing at Q = 27rn, and diverg- 
ing linearly with system size when entering the TCI phase; 
and we consider as well the scaling of the superfluid density 
for both species, Ps,a{b)^ which scales to zero in the thermo- 
dynamic limit when in the TCI phase. Fig. [TT] shows that 



the scaling of both quantities give a consistent estimate of the 
critical jc value for the onset of crystalline order 

The transition to the TCI phase was akeady estimated by 
one of the authors in Ref. [18| for the same density. The es- 
timate of the transition line we present here is corrected with 
respect to that given in Ref. [18J. Indeed in Ref. 1 18 1 the esti- 
mate of jc was biased by the choice of a scaling ansatz on the 
superfluid stiffness which, albeit satisfactorily verified within 
the statistical uncertainties, turns out to not be fully justified. 
The current estimate is not based on any scaling ansatz, and 
it delivers a transition line which is shifted to larger jc values 
(typically by about 10-20 %) with respect to the estimates in 
Ref. ifTSl . 

Moreover Ref. ifTSl identified a narrow region of parameters 
close to the transition line as a super-Tonks regime, character- 
ized by a sub-linearly diverging peak in the density structure 
factor of both species. The subtle effect of the opening of 
a spin gap before entering the crystalline phase was not ob- 
served in Ref. [181, and this led to the erroneous claim of al- 
gebraic off-diagonal correlations for both species in the super- 
Tonks regime. Our current results point to the fact that the 
above regime is rather characterized by superfluidity of both 
atomic species (which persists up to the transition, as shown 



in Fig. 11 1 but with exponentially decaying off-diagonal cor- 
relations, due to the presence of a spin gap. 



